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Abstract 

A low energy non-unitary leptonic mixing matrix is a generic feature of many extensions of the 
Standard Model. In such a case, the task of future precision neutrino oscillation experiments is 
more ambitious than measuring the three mixing angles and the leptonic (Dirac) CP-phase, i.e., 
the accessible parameters of a unitary leptonic mixing matrix. A non-unitary mixing matrix has 13 
parameters that affect neutrino oscillations, out of which four are CP-violating. In the scheme of 
Minimal Unitarity Violation (MUV) we analyse the potential of a Neutrino Factory for determining 
or constraining the parameters of the non-unitary leptonic mixing matrix, thereby testing the origin 
of CP-violation in the lepton sector. 
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I. INTRODUCTION 



There are several indications from particle physics, as well as from cosmology, for the 
existence of physics beyond the Standard Model (SM). For example, the gauge hierarchy 
problem suggests that new physics exists at energies close to the electroweak scale in order 
to stabilise it against large quantum corrections. In cosmology, the evidence for dark matter 
in the Universe requires the extension of the SM particle content. Last, but not least, the 
discovery that neutrinos are massive provides the first clear particle physics evidence that 
the SM has to be extended. 

In general, extensions of the SM will also affect the physics relevant at neutrino oscillation 
experiments. New physics effects on neutrino oscillations are particularly relevant for the 
next generation of precision neutrino oscillation facilities such as Neutrino Factories 
which aim at measuring the unknown leptonic mixing angle ^13, the neutrino mass hierarchy 
[i.e., sgn(Am3]^)), as well as the Dirac phase 6, which can induce CP-violation in neutrino 
oscillations. In most phenomenological studies regarding the sensitivities of future neutrino 
oscillation facilities, the leptonic mixing matrix is assumed to be unitary. 

In contrast to this common practice, it is well known that one generic feature of new 
physics in the lepton sector is the non-unitarity of the low energy leptonic mixing matrix. 
This non-unitarity appears whenever additional heavy particles mix with the light neutrinos 



or their charged lepton partners fs]. After integrating the heavy states out of the theory, 
the 3x3 submatrix of the light neutrinos remains as an effective mixing matrix. This low 
energy leptonic mixing matrix is, in general, not unitary. 

While there are many models of physics beyond the SM which induce non-unitarity, an ex- 
tension of the SM featuring a non-unitary leptonic mixing can be described in a minimal way 
through an effective theory, the so-called Minimal Unitarity Violation (MUV) scheme |4|. 
It contains the relevant low-energy information for neutrino oscillation experiments and is 
minimal in the sense that only three light neutrinos are considered and that new physics 
is introduced in the neutrino sector only. It provides an effective description of all models 
where additional heavy singlets mix with three light neutrinos.^ 



^ Other possibilities to introduce non-unitary leptonic mixing are, e.g., via an additional vector-like lepton 
generation or via fermionic SU(2)l triplets, which are beyond MUV. Non-unitarity in these schemes turns 
out to be significantly more constrained by non-oscillation experiments than in MUV (see, e.g., Ref. Q). 



2 



In MUV, the charged- and neutral-current interactions of the neutrinos {i.e., their cou- 
phngs to the W and Z bosons) are modified. The non-unitary leptonic mixing matrix N, 
which appears in the charged- current interaction, contains the only additional degrees of 
freedom, since the neutral-current interaction of the neutrinos is proportional to A^''"A^ while 
the neutral-current interaction of the charged leptons is unchanged. Thus, instead of the 
three mixing angles and three CP-phases of a unitary leptonic mixing matrix (with only one 
affecting neutrino oscillations), the non- unitary mixing matrix N contains 15 parameters, 
out of which six are CP-violating phases (including two Majorana phases, which do not 
affect neutrino oscillations). 

In this study, wc investigate the potential of a Neutrino Factory for determining or con- 
straining the parameters of the non-unitary leptonic mixing matrix, thereby testing the 
origin of CP- violation in the lepton sector. 

II. GENERAL INTRODUCTION TO UNITARITY VIOLATION 

As motivated in the introduction, non-unitarity of the leptonic mixing matrix is a generic 
manifestation of new physics in the lepton sector. The MUV scheme provides an effective 
field theory extension of the SM and is minimal in the sense that only three light neutrinos 
are considered and that new physics is only introduced in the neutrino sector. Notice that 
this assumption is conservative, since new physics affecting other sectors, such as that of 
the charged leptons, will lead to stronger signals than the ones discussed here. The MUV 
scheme thus describes the relevant effects on neutrino oscillations in the various types of 
models where the SM is extended by heavy singlet fermions (where "heavy" refers to large 
masses compared to the energies of the neutrino oscillation experiments) which mix with 
the light neutrinos. 

In the MUV scheme, the Lagrange density of the SM is extended by two effective 
operators, one of mass dimension five and one of mass dimension six. The dimen- 
sion five operator is the ubiquitous lepton number violating Weinberg operator SC^^^ — 
I c^l^ (^L'^a4>*^ (4>^ L/s^ + H.c, the lowest dimensional effective operator for generating neu- 
trino masses using the field content of the SM. The coefficient matrix cf^^ is of 0{1/M) 
and related to the low energy neutrino mass matrix by rrij, = v^^^c'^^^, where ^ew is the 
vacuum expectation value of the SM Higgs field 0, which breaks the electroweak symmetry. 
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and </) = iT24>*. The SM neutrinos are contained in the lepton doublets La, with a = e, fi,T 
running over the three famihes. 

The effective dimension six operator c^^^ (^La(pj i0 conserves lepton number^ 

and, after electroweak symmetry breaking, contributes to the kinetic terms of the neutrinos. 
After their canonical normalisation, they generate a non-unitary leptonic mixing matrix N, 
as well as non-universal couplings of the neutrinos to the Z boson proportional to N'^N. 
The modified part of the Lagrange density in MUV is given by 

C'^ = ^ {uii^Ui- 17^imiVi+ H.c.) - -1= (1^+ 7^ (1 - 75) A^,, z/^ + H.c.) 

{Z^vn^{l-^^){N^N),^v^ + H.c.) . (1) 
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2 cos 9w 

We note that the MUV scheme is also minimal in the sense that all new physics effects depend 
on the non-unitary leptonic mixing matrix A^. Regarding neutrino oscillation experiments, 
the non-unitarity of A^ affects the processes at the source and the detector as well as neutrino 
propagation in matter. 

To parametrise A^, we use the fact that a general matrix can be written as the product 
of a Hermitian matrix times a unitary matrix. Decomposing the Hermitian matrix as 1 -|- e 
(with e = e^) and denoting the unitary matrix by U , we can write 6| 

N={l + e)U. (2) 

For the complex off-diagonal elements of the matrix e, we use the notation = lea/jje*'^"''. 
Notice that, due to the Hermiticity of e, leapl = \£f3a\ and (pais = —(pfSa- The diagonal 
elements are real and no further parametrisation is required. Constraints on the Eajj can 
also be derived from the experimental data on electroweak decays 0,S|- The present 90 % CL 
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0. 



< 8.0-10-^ l^^^l < 5.1-10-^ 14 and \eee\ < 2.0-10-^ < 



8.0-10"^, I^T-T-I < 2.7- 10""^ [9I]. In our analysis, we will consider unitarity violation consistent 
with the present bounds. Analytic expressions for the neutrino oscillation probabilities in 
terms of U and e can be found in App. \M 

Finally, we would like to comment on other possible parametrisations of a non-unitary 



leptonic mixing matrix. In Refs. 



[Id, 



iil,Q,Q, 



a different parametrisation is advocated, in 



^ We note that since the dimension six operator conserves lepton number, it is not necessarily suppressed 
by the smallness of the neutrino masses. 
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which the deviations from unitarity of the mixing matrix involving the three hght neutrinos 
is related to the mixing between these light neutrinos and the heavy singlets in seesaw type 
theories. The mixing matrix in a seesaw scenario is the unitary matrix that diagonalises the 
extended neutrino mass matrix: 





rriD 1 


f/6x6 = 1 


' m 












Mm j 




\ 


M 1 



(3) 



where rriD and are the neutrino's Dirac and Majorana mass matrices, respectively. In 
the case of only one neutrino family, the unitary matrix is just a rotation of angle 9 ~ mo/M. 
The extension to three or more families is straightforward, performing the diagonalisation 
in two steps: first a block-diagonalisation and then two unitary rotations to diagonalize the 
mass matrices of the light and heavy neutrinos, i.e., 




(4) 



where U and V are unitary matrices. Without loss of generality, we can choose a basis for the 
heavy singlets such that V = I. Analogously to the one family example, when performing 
the block diagonalisation, the mixing between the light and heavy neutrinos is suppressed 
so that 

5 ~ e = mDM^\ (5) 

This suppression is exploited in Refs. jl^ . , where the block diagonalising matrix is writ- 
ten as the product of the 9 possible rotations mixing the light and heavy states and then 
expanded up to second order in the small mixing angles. This results in 

/ 1 I „2 I „2 ^ n n \ 



A 



1 



B 





14'='24 ~T '='15'='25 ■='16''26 2 ^''24 + "^25 ''26/' 
y ■Si4'534 + ^isSgg + S;^gS3g S24"^34 + '^25'^35 + '^26'^36 2 ^'^34 






+ sis + s 



36) 



■^14 -^IS '^16 



'24 ■'25 ■'26 



(6) 



y -534 -535 ■^36 j 

where Sij = Sij exp{i6ij) and Sij 



sin{6ij). Notice that the mixing matrix of the three light 



neutrinos is given by = AU. Thus, the deviation from unitarity, encoded in A, is directly 
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related to the mixing B between the heavy and hght neutrinos. We argue that this is also 
the case with the Hermitian unitarity deviation adopted in Eq. ([2]). Indeed, we can exploit 
the suppression of Eq. to write the unitary block diagonalisation as the exponential 
expansion of an anti-Hermitian matrix: 

I A b\ [ q\ /i-ieet \ 

Thus, the Hermitian deviation from unitarity defined in Eq. is just e = — G0^/2 and its 
relation to the mixing between light and heavy neutrinos in a seesaw scenario is straight- 
forward.^ Furthermore, notice that the deviation from unitary mixing parametrised as in 
Eq. ([6]) can only be applied to the specific case of the mixing between three light and three 
heavy neutrinos while the product of an Hermitian and a unitary matrix is a completely 
general matrix and thus suitable to take into account more general scenarios. In addition, 
the unitarity deviation e is given by the coefficient of the d = 6 operator (e = — 
that modifies the neutrino kinetic terms, introduced in the MUV scheme and obtained in 
the effectwe theory of the seesaw mechamsm after mtegratmg out the heavy smglets (see, 
e.g., Ref. M). 



III. NUMERICAL SIMULATION AND RESULTS 



We will now discuss the sensitivity of future neutrino oscillation experiments to the 
different parameters of the MUV scheme. In particular, we study the Neutrino Factory 



setup proposed in the International Design Study (IDS) [15|, |16J, which consists of and 
z/^ beams from 5 ■ 10^° muon decays per year per baseline. We consider a setting where the 
experiment is assumed to run for five years in each polarity. The parent muons are assumed 
to have an energy of 25 GeV. The beams are detected at two far sites, the first located at 



4000 km with a 50 kton Magnetised Iron Neutrino Detector (MIND) Ii7| and a 10 kton 

n □ 

Emulsion Cloud Chamber (ECC) for r detection ^8|, Il9| , and the second located close to 



the magic baseline 20|, l2l| at 7500 km with an iron detector identical to the one at 4000 km. 



A clean signal of a non-unitary mixing is the presence of "zero-distance effects" stemming 



^ The anti-Hermitian part can be reabsorbed in the unitary rotation, and is thus related to using different 
parametrisations. 
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from the non-orthogonality of the flavour states. Indeed, if the flavour basis is not orthog- 
onal, a neutrino of flavour a can be detected with flavour (5 without the need of flavour 
conversion in the propagation. This translates to a baseline-independent term in the oscil- 
lation probabilities, which is best probed at short distances, since the flux is larger and it 
cannot be hidden by the standard oscillations. For short baselines, this term is [a 7^ /?) 

P,^(L = 0) =4|£„^|2 + 0(£3). (8) 

The oscillation probabilities for longer baselines up to second order in the small parameters 
are derived in App. |Al Near detectors are thus excellent for probing the zero-distance effect, 
in particular r detectors are of importance, since the present bounds on and are 
rather strong. We will therefore study the impact of near r detectors of different sizes 
located at 1 km from the beam source. In particular, we will present all the results for 
near detector sizes of 100 ton, 1 kton, and 10 kton, as well as the results without any 
near r detector. Notice that 10 kton is the detector mass discussed for the ECC detector 
located at 4000 km. However, we have seen no improvement adding such a detector at 
that baseline while the gain in sensitivity that a near detector capable of r detection can 
provide is significant, as we will discuss below. Therefore, we also considered the larger 
mass to show what could be achieved with the planned 10 kton detector located at 1 km 
instead of 4000 km. To simulate the near detector, we use the point-source and far-distance 



approximations. These assumptions are reasonable, 



though somewhat optimistic in the 



high-energy region, as can be seen in Fig. 12 of Ref. However, the loss of flux at higher 
energies, which corresponds to the on-axis neutrinos, may be recovered by using rather 
elongated geometries of the near detector. These are precisely the kind of geometries that 
are being discussed for a magnetized version of the ECC (MECC). Such a detector would 
be limited in size by the above mentioned geometrical considerations and is not likely to be 
larger than 4 kton. On the other hand, all the decay channels of the r could be studied in 
the magnetized version, which would translate into an increase of the efficiency by a factor 
5 with respect to the ECC search for r decays into /i considered here. The impact of near 
/X detectors is still essentially to normalise the neutrino flux and cross-sections, since the 
bounds on and e^e from the unitarity of the CKM matrix and /i — > 67 are particularly 



strong 



In our simulations, we will study the "golden" 



23| and disappearance channels 



in the MIND detectors and the "silver" 



18 



19| i/g — ^ and "discovery" [24^ z/^ chan- 



nels at the ECC detectors, both near and far. For the detector efficiencies and backgrounds, 
we follow the study in Ref. of the MIND detector exposed to the Neutrino Factory beam. 
The efficiencies and backgrounds for the silver channel with an ECC detector are carefully 
discussed in Ref. [3] and we follow the results of that reference. Lacking an analogous study 
for the discovery channel, we assume the same efficiencies and backgrounds as the ones for 
the silver channel described in Ref. 

For our numerical simulations, we scan the complete MUV parameter space, adding nine 
unitarity violating parameters to the six standard neutrino oscillation parameters. The scan 
is performed using the MonteCUBES software 25|, l26||, which allows to perform Markov 
Chain Monte Carlo (MCMC) simulations with GLoBES 27|, |28|]. For the implementation 
of the unitarity deviations in the neutrino oscillation probabilities, we use the NonUnitarity 
Engine (NUE) distributed along with the MonteCUBES package. Using the MCMC tech- 
nique allows the study of possible parameter correlations in the full parameter space without 
restricting the search to varying only a small subset of the parameters. This is due to the 
fact that the number of evaluations required by Monte Carlo techniques increases at most 
polynomially with the number of parameters, while a scan based on grids in the parameter 
space would require to evaluate the event rates and likelihoods at a number of points that 
grows exponentially. For all of our figures, we have used simulations with four MCMC chains 
containing 2 x 10^ samples each. In addition, we have checked that the chains have reached 
proper convergence, in all cases better than R — 1 = 10~^ 2^. It is also important to note 
that, unlike in the standard usage of the GLoBES software, the use of MCMC techniques 
is based on Bayesian rather than frequentist parameter estimation and, as such, the result 
depends on the adopted priors. As priors, we will consider the current bounds on both 
the standard and the unitarity violating parameters, except for parameters to which the 
Neutrino Factory has superior sensitivity, for which we use flat priors. 

Before discussing the more detailed studies, let us comment on some of the general results 
from the simulations. First of all, one of the most remarkable features is that the results do 
not contain significant correlations between any of the unitarity violating parameters, nor 
are the unitarity violating parameters significantly correlated with the standard neutrino 
oscillation parameters. The only exception are some mild correlations between 6*13, 5 and 
the modulus and phase of Ere in the absence of near r detectors which, however, do not 
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lead to new degeneracies between these parameters or spoil the determination of ^13 and S 
at the Neutrino Factory. Furthermore, the addition of a near r detector of only 100 ton is 
enough to almost completely erase these correlations. This implies that the Neutrino Factory 
setup considered here has enough sensitivity to distinguish the effects induced by unitarity 
violation from changes in the standard parameters. Second, the sensitivities of the Neutrino 
Factory to the diagonal parameters of the e matrix, as well as to e^e, do not improve with 
respect to the bounds derived from electroweak decays, which are too stringent to allow for 
observable effects at the Neutrino Factory. Notice that none of the oscillation probabilities 
studied here depend on See, as shown in App. |Al 

We will thus concentrate on the sensitivities to and Ere in the next subsections, 
even though the other unitarity violating parameters and standard oscillation parameters 
are allowed to vary in the simulations. As an example of the sensitivities and correlations 
to all the 15 parameters considered, the 105 projections to the different two-dimensional 
subspaces and the marginalized regions for the 15 parameters can be studied in a triangle 
plot at Ref. {26 1 for the case of no near r detector. The input values chosen for the unknown 
parameters in this example were ^13 = 5°, 6 = 0, \ere\ = 0.005 and (pre = 7r/4, the input 
br the rest of the non-unitary parameters was set to zero. In all our simulations we assume 



30, lail 612 = 33°, 023 = 45°, Aml^ = 8 ■ 10"^ eV^ and Am|i = 2.6 ■ lO'^ eV^. We also 



assumed 4 % priors on 612 and Am^i at la, flat priors were used for the rest of the standard 
oscillation parameters. For the unitarity violating parameters, we consider Gaussian priors 
given by the ranges mentioned in Sec. [Ill 



A. Sensitivity to e-rfj. 

In the left panel of Fig.[Tl we show the sensitivity to the parameter for the four different 
sizes considered for the near ECC. The input values for all the non-unitarity parameters and 
6'i3 were set to zero to derive these curves. We have checked that the results do not depend 
strongly on this assumption. The most remarkable feature of this figure is the extreme 
sensitivity to the real part of which is present already without any near detector. This 
sensitivity mainly originates from the matter effect on the disappearance channel, where the 
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FIG. 1: The 90 % confidence level sensitivity of the IDS Neutrino Factory to the unitarity violating 
parameters fl^ft) o-nd Ere (right). The different curves correspond to different sizes of the near 
T detector, from left to right, 10 kton, 1 kton, 100 ton, no near detector. 

leading non-unitarity correction to the oscillation probability is given by 

= P'Jf - 2 Me,r)AL sin (^^) + 0{e,,), (9) 

where A = \/2GFne, the terms we have omitted here can be found in App. |Al Notice that 
the discovery channel also depends linearly on and that the dependence is CP-violating. 
On the other hand, the mass and efficiency of the ECC detector are much smaller compared 
to those of the MIND detectors for the z/^ disappearance channel and therefore the sensitivity 
is dominated by the latter. As can be seen in the figure, a near r detector will determine the 
modulus of Efj^r through the zero-distance effect. This would translate into a vertical band in 
the left panel of Fig. [T] and thus the increase of the mass of the near detector improves the 
measurement of the imaginary part. However, given the linear dependence due to the matter 
effects on propagation, the bound on the real part from the disappearance channel remains 
stronger. We can also see that the bound on the modulus does not require a very large near 
detector, the bound on the imaginary part is essentially only improved by approximately 
30 % in moving from a 1 kton to a 10 kton ECC detector. 

Another important question is how well the Neutrino Factory would be able to measure 
the unitarity violating parameters if they are non-zero. For this reason, in Fig. O we show 
the sensitivity to assuming that \er^\ = 3.2 • lO^'^ as well as = 45° and —90°, 
respectively, which is disfavoured at only la by current bounds. Thus, this gives a flavour of 
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FIG. 2: The sensitivity of the IDS Neutrino Factory to the unitarity violating parameter Erfi, 
assuming that it takes the value = 3.2 • 10^'^ exp(i7r/4) (left) and e,-^ = — i3.2 • 10^^ (right). 
The different curves correspond to different sizes of the near r detector, from inner to outer curves, 
10 kton, 1 kton, 100 ton, no near detector. 

the best possible situation for actually discovering unitarity violation and a new source of CP- 
violation. Again, we can see that the sensitivity without the near detector is only to the real 
part of Er^. In this setting, there is a degeneracy extending essentially as le^/xl oc 1/ cos(0t-/x), 
along which the real part of is constant and the imaginary part is changing. For the case 
with purely imaginary e^-^ in the right panel of Fig. [21 it is also no surprise that the results 
without the near detector are compatible with = 0. The introduction of near detectors 
results in an effective measurement of I^T^il, i-^-, a vertical band in the plot, which intersects 
the far detector measurement giving rise to two degenerate solutions, one for positive and 
one for negative imaginary part. Again, the actual size of the near detector is not crucial 
and no significant gain is seen beyond 1 kton. 

These figures also show the strong complementarity between the near and far detectors 
when it comes to measuring the phase of the unitarity violating parameter, and thus also 
a non-standard source of CP-violation. Neither the near nor the far detectors alone can 
establish a CP-violating phase by themselves. However, combining the two results excludes 
CP-conservation at 90 % confidence level. 

Note that the slight widening of the allowed region when including the near detector 
results from the use of Bayesian statistics. Since the near detectors discard a large range 
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of allowed values for 0,-^ when \er^\ is close to zero, a slightly larger region in close 
to the correct absolute value of is needed in order to include 90 % of the probability 
distribution. 



B. Sensitivity to 



The right panel of Fig. [T] shows the sensitivity to the unitarity violation parameter Ere 
when the input values for 6'i3 and all the unitarity violating parameters are set to zero. 
Analogously to the sensitivity to Er^, the setup with only the far detectors is more sensitive 
to the real part of the parameter, although the difference is not as pronounced. Furthermore, 
as can be seen in the oscillation probabilities in App. [Aj the probabilities that depend on 
Ere are only the golden, silver and discovery channels, where the dependence is quadratic 
rather than linear, which translates into a weaker bound. Thus, the inclusion of the near 
r detector has a major impact also on the bound which is placed on the real part of Ere- 
Indeed, for a 1 kton near r detector, the sensitivity is essentially flat as a function of (pre 
and is dominated by the near detector. 

Again, the larger mass and efficiency of the MIND detector compared to the ECC trans- 
lates into the golden rather than the silver or the discovery channels dominating the sensi- 
tivity to Ere from the far detectors alone. However, unlike the disappearance channel, the 
golden channel is strongly dependent on the unknown parameters ^^13 and 6 and the input 



values assumed for them will influence the expected sensitivity to Ere- Indeed, the z/g 
probability in presence of non-unitarity is modified to: 

E3L 



efi 



n^,^ + K.psin^ 



+ Im ^Eer 
+ Re <! Eer 



'IE2 , ^ , E^sise'^ 



AL 

sm I I sm 
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E.-A 



-L 



sm 



E.L 



+0 {e') . 



(10) 



where Ei = Amfi/{2E). It is then clear that the relative importance of the real and 
imaginary parts of Ere in this probability strongly depends on the actual values of 613 and 
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FIG. 3: The 90 % confidence level sensitivity of the IDS Neutrino Factory to the unitarity violating 
parameter e-re with 613 = 5° as well as 6 = it /A (left) and 5 = (right). The different curves 
correspond to different sizes of the near r detector, from left to right, 10 kton, 1 kton, 100 ton, no 
near detector. 

5. As an example of this dependence, in Fig. [3l we again show the sensitivity to Sre^ but for 
input values of 6*13 = 5° as well as for 5 = 7r/4 (left panel) and 5 = (right panel). Notice 
that while for 5 = 7r/4 the far MIND detectors are more sensitive to the imaginary part of 
Ere the situation is reversed for 5 = 0. However, the addition of the near r detector for the 
silver channel dominates the bound and the curves incorporating the near detectors forecast 
the same sensitivity regardless of the true values of ^13 and 5. 

In Fig. m we show the analogue of Fig. [2]for ere- In this case, we assume \eTe\ = 5.0 ■ 10~^ 
and (pre = 45° and —90°, which again corresponds to the la disfavoured region. For this 
example, CP-violation would not be discovered for the (pre = 45° case (left panel) at the 
90 % CL, but it would be constrained around its true value already by the far detectors. In 
addition, the inclusion of a near r detector would again constrain the modulus and therefore 
be complementary to the far detector result. For the (pre = —90° case (right panel), the 
complementarity of the near and far detectors is able to exclude CP-conservation at the 
90 % CL. 
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FIG. 4: The sensitivity of the IDS Neutrino Factory to the unitarity violating parameter Ere, 
assuming that it takes the value Ere = 5.0 • 10^'^ exp(i7r/4) (left) and = — i 5.0 • 10~^ (right). 
The different curves correspond to different sizes of the near r detector, from inner to outer curves, 
10 kton, 1 kton, 100 ton, no near detector. 

IV. SUMMARY AND DISCUSSION 

We have considered the sensitivity of the IDS Neutrino Factory setup to minimal uni- 
tarity violation (MUV) by using the Markov Chain Monte Carlo methods implemented in 
MonteCUBES to explore the full parameter space, consisting of the six standard neutrino 
oscillation parameters and nine additional parameters describing the deviation from uni- 
tarity. Our simulations were performed with several different near ECC r detector setups, 
ranging from no near detector to near detector masses up to 10 kton. 

Our results imply that the Neutrino Factory will be excellent for probing some of the 
unitarity violating parameters. In particular, a sensitivity of (9(10~^) to the real part of 
the unitarity violating parameter is found. This is mainly due to the matter effects 
in the disappearance channel at the far detectors, for which the oscillation probability 
is only linearly suppressed in Re(£:^^). On the other hand, we find that a near r detector 
with a mass as small as 100 ton would dominate the sensitivity to Sre-, as well as that to 
the imaginary part of Sr^., through the measurement of the zero-distance effect, providing 
sensitivities down to 0(10"'^). For the other unitarity violating parameters, we recover the 
priors of our simulation, which were set to the current experimental bounds. The setup 
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studied here will therefore not improve our present knowledge of them. 

Furthermore, we find no degeneracies neither among the different unitarity violating pa- 
rameters, nor between the unitarity violating parameters and the small standard neutrino 
oscillation parameters, such as ^13. This means that the sensitivities to the standard oscil- 
lation parameters are robust even in presence of unitarity violation. 

Regarding the prospects of an actual detection of unitarity violation, and especially CP- 
violation stemming from non-unitary mixing, we find that the near and far detectors play a 
very complementary role. In the case of Erfj,, the far detectors are only sensitive to the real 
part of the unitarity violating parameter while the near detector can measure its modulus, 
neither is sensitive to unitarity violating CP-violation by themselves. However, it can be 
effectively probed by considering the combination of the two, as illustrated in Fig. [2J 

ivity to unitarity violation at a Neutrino 
, the sensitivity to the real part of due 
to matter effects has not been discussed (however, a similar term in the z/^ disappearance 
channel is present in and has been studied for the case of oscillations into sterile neutri- 
nos {2^). Furthermore, these studies have not systematically scanned the parameter space 
while keeping all parameters free within their prior values. Thus, the observation that 
there are no extended degeneracies, neither between the standard and unitarity violating 
parameters, nor among the unitarity violating parameters themselves, is also new. 

We conclude that a Neutrino Factory would provide powerful tool for probing unitarity 
violation in the leptonic mixing matrix. For the parameters to which it is most sensitive, the 
sensitivity is an order of magnitude better than the current experimental bounds. Finally, 
the interplay between the near and far detectors would allow to test new sources of CP- 
violation in the lepton sector. 



We would like to stress that, w. 
Facto. Ka. bee,. .ud,ed«o.fl 



lile the sensi 
d, 33, 34, 35 



Acknowledgments 

We would like to thank A. Donini, W. Wang, and Z.-z. Xing for useful discussions. 
This work was supported by the Swedish Research Council (Vetenskapsradet), contract no. 
623-2007-8066 [M.B.]. S.A., M.B. and E.F.M. acknowledge support by the DFG cluster of 
excellence "Origin and Structure of the Universe". J.L.P. acknowledges financial support 
by the Ministry of Science and Innovation of Spain (MICIINN) through an FPU grant. 



15 



ref. AP2005-1185. S.A., E.F.M. and J.L.P. also acknowledge support from the European 
Community under the European Commission Framework Programme 7 Design Study: EU- 
ROnu, Project Number 212372. 



APPENDIX A: OSCILLATION PROBABILITIES IN THE PRESENCE OF UNI- 
TARITY VIOLATION 

In this Appendix, we derive the probabilities Pajs in matter assuming constant density. 
In order to perform the calculation, we will use the Kimura-Takamura-Yokomura (KTY) 
fonngsm flH. which has already bee„ apphed to the MUV scheme the Appendix of 
Ref. [6|. Since the constraint on is strong enough to safely neglect Ee^ in the oscillation 
probabilities, we will not consider it below. However, it has been considered in the numerical 



analysis presented in the main part of this paper. The effective flavour eigenstates are given 
by: 



/aa\ 



The parameters that appear linearly in the normalisation factors are Eee, Sfifi, and Err, which 
are already better constrained by other considerations than the sensitivities we find for a 
Neutrino Factory. Thus, the determination of the fluxes and cross-sections by the near 
detectors only suffer from a minor additional theoretical uncertainty. We will present the 
oscillation probabilities P^Ua ^/s) = Pa/s without taking the normalisation factors into 
account. Notice that this will not be at all relevant for the golden and silver channels, since 
the probabilities are already order before taking the normalization factors into account. 
Thus, the corrections would be at most 0{e^). 

The oscillation probability Pajs, expressed as a function of the KTY parameters, is 

P.P = \{NN^U\^-Aj2MXfxr)sin^ 

j<k \ 

+2 J2 MXfxf*) sm{AEjkL), (A2) 

j<k 

where AEjk = E^-Ek and Xf = {N*W)ajiN*W)lj (j = 1, 2, 3). Here, Ei are the effective 
eigenvalues in matter and Wij is the unitary matrix which diagonalizes the evolution equation 
for the mass eigenstates: 
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SM expansion parameters (rj) 


MUV expansion parameters 


013, AmlJAmj-^, 6623 = 023 - vr/4 





TABLE I: The small expansion parameters used in our neutrino oscillation probabilities. We will 
refer to the set of SM expansion parameters as ij. The full set of expansion parameters will be 
referred to as e, while only the set of MUV expansion parameters will be denoted by Sap- 



(A3) 



where Ei = ^m1^/ {2E). Assuming that the electron and neutron number densities are 
equal^ (z.e., rie = Un), Ti. can be expressed as 



AAA 



n = diag(Ei, E2, Es) + Art ^iag j N, (A4) 

where A = \piGpnf,. Finally, according to the KTY formalism applied to the MUV scheme 
(again, see Ref. [6|), X"^ can be expressed as 



(A5) 



where 



(A6) 



(AE2iA^3i)~' (^2^3, -E2 - Es, 1) ^ 
-{AE2iAEs2)~^ {E3E1, —Es — Ei,l) 
y {AEsiAEs2)-' {E2E,, -E2 - El, 1) ^ 
Once the effective eigenvalues in matter are known, it is straightforward to obtain the 
expressions for the neutrino oscillation probabilities. However, in order to obtain reasonably 
simple expressions, it is necessary to expand them in small parameters. Here, we present 
the oscillation probabilities to second order in the parameters listed in Tab. [B 

To second order in e, we can find the eigenvalues by using perturbation theory. We find 
that 



El = A 



A 



4- f -L _ 

A-E3 2 2 



+ Oie') , (A7) 



This is a very good approximation in the case of neutrino oscillations in the Earth. 
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A 



E2 ^2 
A 
1 , 



^12 4^2 



sin2(2^i2) +Re(£^^ 



1 + + 



1. 
2 



2 



^3 



4t 

4'^ 



4^3 

1 



5^23 
^2 ^ + 



'MM' 



72 



(A8) 



1 + Iji^i^t^ 



fj.fl + ^Tt) ^ [^TT ~^ ^HfJ.) 



- flT I 



- + 



+ 56*23 (^rr 



/I V rr ' fJifi/ r\ ' A \ , (A9) 

Notice that, for — > 0, we recover the SM results as expected. These resuhs allow us 
to obtain V^^ at second order. Thus, we only need to compute Y^"^ at the same order, 
the computation is straightforward but tedious (see Eq. (lASP ). For brevity, we do not 
present the results for V'"^ and Y^^ here. However, we would like to comment that, for 
the golden and silver channels, it is enough to compute these quantities to first order, since 
is already of first order in 77. This is not true in the case of the v,-Vr sector, where 



-^3 le=0 



|£ = 



X. 



Tfl\ 



k=0 



1/2. The advantage of this sector, from the 



point of view of discovering new physics, is that the effects of the new physics can appear in 
the probability at first order as an interference term between the SM and the new physics 
without additional suppression by 77. For this reason, we keep only the interference between 
the O {Safi) terms and the O (77) ones at second order^ in that sector. 

In the end, we obtain the following expanded oscillation probabilities at the orders men- 
tioned above: 



^lfl 



fifi • ' ""MM 



269. 



A 



23[£ 



MM 



sm 



'"E^j \ 2 
[2 Reiser) - msiSfif. - £rr)] AL siniEsL) + O [elp) 



(AlO) 



P, 



fXT 



flT 



MT| 



+ 



2 Re(£:^^ + Err) + 856'23(£:^^ - £ 



A 



sm 



E^L 



^ It could also be justified to neglect the O I Safj— — ^ ) terms, since the maximal allowed value of — 



is at least one order of magnitude smaller than the maximal allowed values of S13 and (5023- However, we 
keep also these terms for completeness. 
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+ [-2 lm{ef,r) - 5e2^{e„„ - erT)AL\ siniEsL) 



-v^Imjee 
+V2Re\ee 



E2 . . . (AL\ fEs-A 
— sin(2^i2) sm ( — 1 cos ( — - — L 

2^35136'^ fAL\ . fE^-A^ 

■ — — cos sm L 

A-E. \ 2 \ 2 



AL\ . /E3L\ . fEs-A^ 
sm I — — sm — — sm L 



+0 (elfs) , 



sm 



(All) 



Peu — Ppii ~l~ l^erl Sm 



e/j. 



en 



E.L 



+ Re <^ £er 



--^sm(2^.) + -^-^ 



AL\ 
sm I sm 

2 ; 



+0 ie') , 



1 E2 . ^ . AL\ E,- 
--sm(2^,,)sm(^— jcos^^ 

2^/2 E: 



-^—^ cos ( — Ism f^-L 



sm 



E^-A 



sm 



(A12) 



pSM , ^1 



-2 
— Im 
-2v^Re<!e 



2-2 



2 ^^^^^^ Re(£e.e'^ 



sm(2^i2) Re(£eT) 



sm 



2 / ^3 - ^ 



sm 



-iS 



sin(2^,2) sin(AL) - ^^f^^" sin({£;3 - 



1^2 . X -^35136 

2T""('^-)-:4^ 



+ Im l^e^ sin(26'i2) sin ^"^"^ cos ^ 



sm I I cos I I sm I L 



+ 



2^/2E: 



35136"* fAL\ . [Es-A 

— — cos sm L 

A-E, \2 \ 2 



cos 



E:>L 



+0 {e') 



(A13) 



Notice that we do not neglect the zero-distance effect in the v^j^-Vt sector. Although this is 
not within the order of the expansion, we keep it as it plays an important role in the analysis 
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of the neutrino flavour transitions at near detectors. 
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